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Abstract 

We investigate the quark mass dependence of the nucleon mass M^. An interpo- 
lation of this observable, between a selected set of fully dynamical two-flavor lattice 
QCD data and its physical value, is studied using relativistic baryon chiral pertur- 
bation theory up to order p 4 . In order to minimize uncertainties due to lattice dis- 
cretization and finite volume effects our numerical analysis takes into account only 
simulations performed with lattice spacings a < 0.15 fm and m n L > 5. We have also 
restricted ourselves to data with m n < 600 MeV and m sea = m va \. A good interpo- 
lation function is found already at one-loop level and chiral order p 5 . We show that 
the next-to-leading one-loop corrections are small. Prom the p 4 numerical analysis we 
deduce the nucleon mass in the chiral limit, Mq ~ 0.88 GeV, and the pion-nucleon 
sigma term utv = (49 ± 3) MeV at the physical value of the pion mass. 



1 Work supported in part by BMBF and DFG. 



1 Introduction and framework 



Lattice QCD on one side and chiral effective field theory, on the other, are progressively 
developing as important tools to deal with the non-perturbative nature of low-energy QCD 
and the structure of hadrons £Q . The merger of both strategies has recently been applied 
to extract physical properties of hadrons — such as the nucleon — from lattice QCD simula- 
tions. Of particular interest in such extrapolations is the detailed quark mass dependence 
of nucleon properties. Examples of recent extrapolation studies concern the nucleon mass 
0E!, its axial vector coupling constant and magnetic moments [HE], form factors [Hj and 
moments of structure functions [7] . 

Accurate computations of the nucleon mass with dynamical fermions and two active 
flavors are now possible [3111021 i n lattice QCD. However, the masses of u— and d— quarks 
used in these evaluations exceed their commonly accepted small physical values, typically 
by an order of magnitude. It is at this point where chiral effective field theory methods 
are useful - within limitations discussed extensively in refs. (21 IHJ - in order to interpolate 
between lattice results, actual observables and the chiral limit (m u ^ — >• 0). In this paper 
we explore the capability of such an approach for extracting the nucleon mass and the 
pion-nucleon sigma term. 

The nucleon mass is determined by the expectation value (N |0^| N) of the trace of the 
QCD energy-momentum tensor [TT] . 

% = ^r-G^G*™ + m u uu + m d dd+ ... , (1) 

where G^ u is the gluonic field strength tensor, (3(g) is the beta function of QCD and m q qq 
with q = u, d . . . are the quark mass terms (we omit here the anomalous dimension of the 
mass operator, as in [TJj ). So the physical nucleon mass M N can be expressed as 

M N = Mq + ctn (2) 
in terms of its value Mo in the SU(2)f chiral limit, 

M = (N\^-G, V G^ + ...\N) (3) 

(with suitably normalized nucleon Dirac spinors). The dots refer to possible contributions 
from heavier quarks, other than u and d, and the sigma term is defined as 

a N = m q—, — — = (N \m u uu + m d dd\ N) . (4) 
/ — ' dm„ 



q=u,d 



'i-l 



The quark mass dependence of translates into a dependence on the pion mass: ~ m q 
at leading order. We pursue this connection in the symmetry breaking part of the chiral 
effective Lagrangian. 
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The framework of our study is relativistic SU(2)f baryon chiral perturbation theory 
(BChPT) as described in ref. [HI]- The effective Lagrangian required for our analysis of the 
nucleon mass up to 0(p 4 ) is 



C = £g>+£g>+£« + 4 2 > (5) 

with 

= Cl Tr(x + )M-^Tr(^)(W^ + h.c.) + ^Tr(«X) + 
6§ = e 38 (Tr( X+ )) 2 M + ^Tr( X ^- X 2 _)^ 

(Tr( X 2 _) - (Tr( X _)) 2 + Tr( X 2 + ) - (Tr( X+ )) 2 ) + (6) 

In C$ we follow the notation of ref. 14J. Here is the leading order pion Lagrangian 
including the mass term. The nucleon Dirac field is denoted by and M is the nucleon 
mass in the chiral limit. The axial field w M and the covariant derivative _D M involve the 
Goldstone boson fields via U(x) G SU(2), and x± — u*x u ^ ± u X^ u i v 2 — U, parametrizes 
the explicit chiral symmetry breaking through the quark masses; here we use x = 2BA4, 
where M, = diag(m u ,m d ) and B = —(qq)/f% is the chiral condensate divided by the pion 
decay constant squared, both taken in the chiral limit. In the following we neglect isospin 
breaking effects. 



2 Analytic results 

2.1 0{p 3 ) Analysis 

The leading order contribution to the shift of the nucleon mass from its value in the chiral 
limit comes from the explicit chiral symmetry breaking piece in C N , which drives the nucleon 
sigma term of Eq.(4). The next-to- leading order (NLO) contribution is represented by 
diagram (a) of FigUl with the irNN vertex generated by C% ■ We have evaluated the relevant 
one-loop integrals using the so-called infrared regularization method jT^j. It represents 
a variant of dimensional regularization which treats one-loop integrals involving baryon 
propagators in a way consistent with chiral power-counting. The diagram (a) develops a 
divergence proportional to m\. It is absorbed in contact terms which are formally of fourth 
order. We denote 2 the counterterm structure that renders the 0(p 3 ) contribution finite 
as — eim^^S . In the notation of ref. ^3] it involves the coupling constant combination 
e 1 = -(16e 38 + 2&L15 + 2e U6 ) from dff . 

Following the reasoning outlined here, the constraint to obtain a finite result at leading 
one-loop order has effectively promoted a linear combination of p A couplings — denoted by 

2 Our coupling ei differs from the convention of ref.p] by a factor 4. 
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Figure 1: One-loop graphs of NLO (a) and NNLO (b, c) contributing to the nucleon mass 
shift. The solid dot denotes a vertex from , the diamond a vertex from C$ . 



e\ — into the p 3 calculation. The resulting expression for the m n dependence of Mat then 
reads: 



M 



N 



M - 4 Cl ml + 



eJ(A) + 



64vr 2 / 2 M 



2hA 
A 



3<M 2 



4M 2 



7r rat 
— h arctan — - 
2 \/4M 2 m2 - mi 



167T 2 / 2 ^ 

Here e^(A) is the finite (renormalization scale A dependent) part of e\ 

3 9A 2 T 



(7) 



ei 



eJ(A) 



and any ultraviolet divergences appearing in the limit d — > 4 are subsumed in 



A 



d-i 



167T 2 



^ I -^(ln(4vr)+r'(l) + i; 



For further discussion we expand the 0(p 3 ) result Eq.((7|) in powers of the pion mass and 
obtain 



M, 



N 



M - 4 Cl ml 



:ra. 



+ 
+ 



e[(X)- 



3^ 2 



3^ 2 

327T/2"" 

3^ 2 



64^Af/ 1 + 2hl -» 



2567r/ 2 M 2 w 



0(m«). 



Note that the sum of the first three terms in this formula coincides with the well-known lead- 
ing one-loop expression for of Heavy Baryon Chiral Perturbation Theory (HBChPT), 
as expected in the infrared regularization approach From Eq.fjHJ) one can also deduce 
that the counterterm e\ of Eq.fJJJ), required in relativistic baryon ChPT for renormaliza- 
tion purposes, is equivalent to the counterterm introduced in j3] which regularizes the short 
distance behaviour in HBChPT. 



In ref.jH] an assessment of convergence properties of ChPT has been performed. It was 
shown that, up to pion masses around ~ 300 MeV, the chiral perturbation expansion de- 
velops a stable plateau region independent of cutoff scales. Moreover, this analysis indicates 
that an upper limit for this plateau behaviour may be reached when approaches about 
600 MeV. While these considerations were made in a non- relativist ic framework, explicit 
comparison shows that the relativistic approach used in the present work contains the same 
chiral structures as those discussed in ref. [3]. We can therefore assume that, with respect to 
the internal consistency of ChEFT, our analysis is applicable for pion masses well above the 
physical one. We give m n < 600 MeV as an estimate for the range of validity. We emphasize 
that, in contrast to the framework adopted by the Adelaide group j2], all the terms beyond 
the leading c\ contribution to the nucleon mass in Eq.fjHJ) are part of the same chiral order 
p 3 . The numerical evaluation of the individual contributions to Eq.fjHJ) (see Section 3) shows 
that the large fluctuations in the chiral extrapolation using dimensional regularization, re- 
ported in [2], arise from examining only the first four terms in Eq.fjHJ), instead of keeping 
the full expression (JZj). 



2.2 0{p A ) Analysis 



Let us now focus on the next-to-NLO (NNLO) contribution to the pion mass dependence 
of the nucleon mass. This involves also graphs (b) and (c) in FigQ which include vertices 
generated by as well as wave- function renormalization In order to avoid having 

to deal with a number of counterterms too large to be handled in a meaningful numerical 
analysis, we decide to truncate the chirally expanded formula at 0(m^). We will show 
numerically that this truncation approximates the full function reasonably well for the 
parameter ranges considered here. Up to terms of order m% no counterterms other than e\ 
are required for a finite result. At 0(p A ) one then obtains 



M 



N 



M - 4 Cl ml 



+ 



3 9 a 2 

327T/2 

(g_l 

M 



3^ 2 



256tt/2M ( 



5 1 0(m 



e i( A ) - ^7Zm \ — ~ 2" 



32tt 2 / 2 \M 



^- - 8ci + c 2 + 4c 3 ) In ^ 



A 



■m 



(9) 



where now 



ei 



e[(X) 



3L 



9aT_ 
Mo 



lCi + C 2 + 4C; 



This expression includes the constants c 2 and C3 which encode the influence of the A (1232) 
resonance in low energy pion-nucleon scattering. The terms up to m\ have already been 
discussed in [TE\. For related discussions of the SU(3)f case see jTH] and references therein. 
It is also interesting to observe that our truncation of the relativistic result at m\ as shown 
in Eq.(j2J) formally coincides with the expansion of nucleon mass in HBChPT to fifth order, 
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as there are no genuine two-loop graph contributions at this order in the chiral expansion 

H3. 



3 Numerical analysis and contact with lattice QCD 

We proceed with the numerical evaluation of Eqs.(J7J) and 0. We set the nucleon axial 
vector coupling and the pion decay constant equal to their physical values, g& = 1.267 and 
f n = 92.4 MeV. Strictly speaking, these quantities should be taken in the chiral limit. We 
have checked that using current estimates for g° A and /° at m n — > does not lead to any 
significant changes in our final results. Details are discussed in the last part of this section. 
Without loss of generality we choose A = 1 GeV. At order p 3 we are then left with three 
unknown parameters (Mo, c\ and e^(lGeV) = ei) and four parameters at order p A (Mo, 
ci, A = eJ(lGeV) + 3c 2 /(128tt 2 ^) and B = c 2 + 4c 3 ). Our NNLO result is identified with 
Eq.(JHJ). This limits the number of tunable coefficients but still keeps sufficiently many orders 
in to provide a good approximation to the full 0(p 4 ) result. 

The unknown parameters are determined using as input a combined set of lattice QCD 
data obtained by the CP-PACS 0, JLQCD and QCDSF [TU] collaborations. These com- 
putations are performed using fully dynamical quarks with two flavors. In order to minimize 
artifacts from discretization and finite volume effects, we have selected from the whole set 
of available data those with lattice spacings a < 0.15 fm and L > 5. Furthermore we re- 
strict ourselves to the resulting four data points with m n < 600 MeV and with equal valence 
and sea quark masses, m sea = m va i. A study of finite volume dependence is in preparation 
[T9"] . We have expressed lattice data in physical units via the Sommer scale r = 0.5 fm [2Hj . 
not taking into account systematic errors arising from possible quark mass dependence of 
r occuring in dynamical simulations. A forthcoming study will address this issue |2*T] . 

In a preliminary step we have fitted the set of lattice points using the LO result treating 
M and C\ as free parameters. The resulting linear fit gives an estimate of C\ about a factor 
3 smaller than the value determined in vrA^ scattering analyses with Xdof = 2.25. We 
conclude that linear fits in the quark mass are not appropriate to describe the quark mass 
dependence of baryon masses. 

Successive steps in our analysis are shown in FigsHJ HU and summarized in table H We 
have first analyzed the 0(p 3 ) result, Eq.(J7J) (Fit I). The best Fit I curve is the solid one 
drawn in FigJ21 The low-energy constants come out of natural size. Furthermore, C\ which 
determines the slope of M^m 2 ) for small has the correct sign and the value of t\ is 
within the range quoted in ref. [3]. 

As seen in Fig|2J the curve obtained by fitting the four lattice data with m n < 600 MeV 
and including the physical point as a constraint shows a surprisingly good (and not yet 
understood) agreement with lattice data even up to m n w 750 MeV. 

The same figure also shows how Fit I develops term by term when the full order p 3 NLO 
expression (J2J) is expanded according to Eq.(jHJ). We emphasize again that in the hierarchy 
of terms with increasing powers of m n , as represented by the dash-dotted, short-dashed and 
long-dashed curves, all contributions are of the same chiral order p 3 in the formulation of 
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Figure 2: Solid/dotted line: best fit curve using the one-loop result at chiral order p 3 Eq.(J7J). 
Input: four lowest lattice data points with m n < 600 MeV and physical nucleon mass (Fit 
I). The dotted extension of this curve for m 2 > 0.4 GeV 2 indicates the region where the 
application of baryon ChPT is usually believed to become unreliable. For illustration we 
also show a subset of lattice data up to m n rs 0.8 GeV, those compatible with the cuts in 
lattice spacing and volume as explained in the text. The solid dots are CP-PACS data, 
the boxes refer to JLQCD and the empty circles to QCDSF. The dot-dashed, dashed and 
long-dashed curves show, respectively, the contributions from the sum of the first three, four 
and five terms in Eq.©. 



baryon ChPT we use. Evidently, truncating the expansion (jSJ) at m% already provides a 
decent approximation to the full 0(p 3 ) result. 

In the NNLO case the statistics of our restricted data sample is not sufficient to constrain 
all the parameters. We have therefore used input values for C2 and C3 available in the 
literature. We set C2 = 3.2 GeV -1 in agreement with refs. |22[ |2"3~] and performed two kinds of 
fits, one with C3 = —3.4 GeV -1 , found in j23] to be consistent with empirical NN phase shifts 
and still within the error bar quoted in [23] (Fit II), and another one with C3 = —4.7 GeV" 1 , 
the central value determined by [2H] in low-energy nN scattering analysis (Fit III). 

Fit III underestimates c\, whereas the value obtained in Fit II for this LEC is in agree- 
ment with and with the outcome of the analysis by Becher and Leutwyler of low-energy 
irN scattering in SU(2) f relativistic baryon ChPT [2E], the framework we use. Furthermore, 
the value for C3 employed in Fit II is quite close to —2.9 GeV" 1 , the one corresponding to 
the empirical spin-isospin averaged p-wave scattering volume which is dominated by the 
A (1232) contribution. 

FiglSl demonstrates that the difference between the 0(p A ) and 0(p 3 ) results is relatively 
small over the entire range of that we analyzed. We explicitly show that higher order 
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Figure 3: Solid curve: best fit (Fit II) using the NNLO result, Eq.Q, at chiral order p A . 
Dashed curve: NLO result (chiral order p 3 ) from Eq.(J2J) using as parameters the central 
values of Fit II. e± has been deduced setting C2 = 3.2 GeV" 1 . For details on the data points 
see Fig|21 



chiral corrections are small, even at pion masses well above the physical one. We therefore 
believe that our interpolation has passed the necessary tests of consistency and convergence 
for m„ < 0.6 GeV. 

In the calculations underlying the Fits I-III we have used g& = 1.267 and /„. = 92.4 MeV 
as input. Rigorously speaking, we should have used values of those quantities in the chiral 
limit. We have performed test calculations with g\ = 1.2 jS] and /° = 88 MeV ^H]- With a 
slight re-adjustment of e\ by less than 3%, any one of the quantities in table Q changed by 
less than 1% when replacing f n with g Q A , f®. 



3.1 The sigma term of the nucleon 

The pion-nucleon sigma term ctjv, defined in Eq.(j3J, translates into 

2 9M N 

0-jv = m^j—- , 10 

if we assume that the Gell-Mann - Oakes - Renner (GOR) relation ~ m q holds and 
we can neglect 0(m^) terms. An improved analysis of s-wave tttt scattering lengths (2E1 
indicates that the 0{m 2 q ) corrections to the GOR relation are very small, although this 
statement becomes progressively less accurate with increasing quark masses, and further 
detailed examination of the role of strange quarks in this context is necessary. A recent 
systematic analysis of results for pseudo-Goldstone boson masses from Nf = 2 lattice 
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m^IGeV 2 ] m^IGeV 2 ] 

Figure 4: The pion- nucleoli sigma term as a function of from Eq. (jl2jl . using as input 
the central values from Fit II (see table Q). The small m n region is magnified in the right 
panel and plotted together with the frequently quoted empirical = 45 ± 8 MeV |3D] . 



QCD comes to conclusions consistent with those drawn in ref.j2Hj- We therefore use Eq. fjlO)) 
in the following. 

At chiral order p 3 starting from Eq. (jHJ one finds the expression 



2 %a 3 



a N = -A Cl m n - —r—^rn^ 
64tt/^ 



2 



The corresponding NNLO result is derived from Eq.(jOJ): 

A 2 99A 2 3 



+ 



4 



W\ ~ ufe (S - 6C1 + 3C3 ) - life (f - 8C1 + C2 + 4C3 ) ln T 

^ H + 0«). (12) 



Our deduced values of ojv at the physical point are summarized in table El The behaviour 
of the sigma term as a function of the pion mass is shown in FigJU Within errors, this curve 
is compatible with the "empirical" sigma term = 45 ±8 MeV extracted in ref.|3U], but it 
does not favor the much larger value reported in ref.jHJ. O ur result is also consistent with 
the analysis of ref.jH2], within the larger uncertainties quoted there. 
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4 Discussion and conclusions 



The present work has been aimed at improving and updating interpolations of the nucleon 
mass, using chiral effective field theory, between the range of relatively large quark masses 
accessible in full lattice QCD simulations, and the small quark masses relevant for compari- 
son with physical observables. A remarkably good interpolation can already be achieved by 
a one-loop calculation at chiral order p 3 using relativistic baryon ChPT. In either case short 
distance dynamics, including effects of the A(1230) and possibly other resonance excitations 
of the nucleon, are encoded in a single counterterm that controls the contributions of order 
m\. 0(p 3 ) relativistic baryon ChPT is therefore free of limitations of HBChPT discussed 
in [21 E]- Apart from the nucleon mass in the chiral limit, the only remaining parameter, Ci, 
drives the pion-nucleon sigma term. Our interpolation is based on a selected set of lattice 
data corresponding to the largest available lattice volumes and the lowest available pion 
masses, in order to minimize uncertainties from finite size effects and from quark masses 
too large to be handled using perturbative chiral expansions. Surprisingly, the resulting 
interpolations work even in a pion mass region where the approach is commonly believed to 
become unreliable. 

The extension to NNLO (chiral order p 4 ), truncated at order m%, introduces in addition 
the pion-nucleon low-energy constants 02,3 which primarily reflect the impact of A resonance 
physics on low-energy irN dynamics. We have constrained the input values for these two 
LECs from irN phenomenology. The fit interpolating the nucleon mass between the chiral 
limit and the lattice data remains remarkably stable and even improves slightly when going 
from NLO to NNLO. Our analysis explicitly shows that 0(p 4 ) corrections are small with 
respect to the 0{p 3 ) result. 

The pion-nucleon sigma term deduced from the dependence of the nucleon mass in 
the NNLO "best fit" (Fit II) is fully consistent with that obtained by Gasser, Leutwyler 
and Sainio jSU]- This is a nontrivial result since no such constraint has been built into the 
procedure. 

In summary, the outcome of the present study is promising. It demonstrates that ex- 
trapolation methods based on chiral effective field theory can be successfully combined with 
lattice QCD results in order to bridge the gap between simulations and observables. Of 
course, remaining uncertainties need to be further investigated, such as corrections due to 
finite lattice volume and questions concerning convergence properties of the chiral expansion 
with quark masses exceeding 100 MeV. 

Future studies will include the quark mass dependence of r , Mjv//tt and the implicit 
quark mass dependence of f n and q-kNN- 



We gratefully acknowledge many stimulating discussions and communications with M. 
Birse, M. Gockeler, H. Leutwyler, U.-G. MeiBner, G. Schierholz and A.W. Thomas. We 
thank the QCDSF-UKQCD Collaboration for providing us with their data prior to publi- 
cation. 
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Table 1: Fit results for MN(m n ) described in detail in the text. Fit I refers to the inter- 
polation based on the 0(p 3 ) NLO result, Eq.(J7J). Fit II and Fit III are based on the 0(p 4 ) 
NNLO result, Eq.©, respectively with c 3 = -3.4 GeV" 1 [21] and c 3 = -4.7 GeV" 1 |2S]. 





M [GeV] 


a [GeV- 1 ] 


ei [GeV -3 ] 


A [GeV~ ;i ] 


B [GeV" 1 ] 


Fit I 


0.891 ±0.004 


-0.79 ±0.05 


3.5 ±0.6 






Fit II 


0.883 ±0.003 


-0.93 ±0.04 




3.8 ±0.6 


-10.4 (fixed) 


Fit III 


0.872 ±0.003 


-1.11 ±0.04 




4.1 ±0.6 


-15.6 (fixed) 



Table 2: The pion-nucleon sigma term deduced from the NLO and NNLO fits for Mn{itl. k ) 
given in Table 1. 





a N [MeV] 


Fit I 


43 ±4 


Fit II 


49 ±3 
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